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Nonasymptotic  Formulae  for  the  Distribution 
of  the  Maximum  of  Smooth  Gaussian 

Processes 


Jean  Diebolt  and  Christian  Posse 


Abstract 

We  derive  an  integral  formula  for  the  density  of  the  maximum  of  smooth 
Gaussian  processes.  This  expression  induces  explicit  lower  and  upper  bounds 
which  are  in  general  asymptotic  to  the  density.  Our  constructive  approach 
relies  on  a  geometric  representation  of  Gaussian  processes  involving  a  unit 
speed  parameterized  curve  embedded  in  the  unit  sphere.  ' 


1  Introduction 


Let  X(i),  t  €  /  =  [0,r]  be  a  real  Gaussian  process  with  mean  zero  and  continuous 
sample  functions.  Numerous  papers  have  been  devoted  to  the  study  of 

Z  =  supX(t) 

tel 

(see,  e.g.,  Adler,  1981,  1990,  Leadbetter  et  a/.,  1983,  Samorodnitsky,  1991,  for 
a  survey).  It  turns  out  that  the  exact  distribution  of  Z  is  known  only  for  the 
Wiener  process,  the  Brownian  Bridge  B(t),  B{t)  -  J5(u)dii  (Darling,  1983),  the 
Ornstein-Uhlenbeck  process,  the  integrated  Wiener  process  (see,  e.g.,  LachaJ,  1991), 
a  class  of  “sawtooth”  processes  (see,  e.g.,  Cressie,  1980)  and  the  random  cosine  wave 
X(t)  =  ^1  cos(t)  -I-  ^2  sia(t),  where  and  ^2  are  independent,  identically  distributed 
A/^(0, 1).  Otherwise,  the  results  obtained  are  either  asymptotic  formulae  for  the  tail 
probabilities  F{Z  >  a}  as  a  00  or  nonasymptotic  upper  bounds  for  P{Z  >  a}. 

*AMS  1991  subject  classifications.  Primary  60G15,  60G70.  Secondary  60G17. 

Key  words  and  phrases.  Differential  geometry,  Gaussian  processes,  extreme  value,  nonasymptotic 
formulae,  density. 
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In  this  paper,  we  consider  the  class  of  centered  Gaussian  processes  of  the  form 

^(0  =  2<n<DC,  (1) 

where  j  >  1,  are  independent,  identically  distributed  A'’(0, 1), 
and  the  functions  Tx{t)  >  0  and  gj{t),  j  >  1,  are  sufficiently  smooth. 

The  existence  of  the  Karhunen-Loeve  expamsion  for  Gaussian  processes  with 
continuous  sample  functions  ensures  that  this  class  is  very  large  (Adler,  1990). 
First,  all  Gaussian  processes  with  smooth  covariance  function  have  such  a  rep¬ 
resentation.  Second,  more  general  processes  can  be  arbitrarily  approximated  by 
processes  of  this  class  with  respect  to  the  uniform  norm.  In  addition,  in  the  context 
of  stationary  Gaussian  processes,  the  discretization  of  the  hexmonic  representation 
^(0  =  exp(itA)dA(A)  (where  A  is  a  Gaussian  measure  with  independent  in¬ 
crements)  or  of  the  filtered  white  noise  representation  X{t)  =  f^h{t  -  s)dlF'(s) 
(where  W{t)  is  the  Wiener  process)  induces  the  form  (1).  See  also  Berman  (1988) 
for  stochastic  modelizations  leading  to  processes  (1). 

Berman  (1988)  studies  the  asymptotic  behavior  of  the  tail  probabilities  of  the 
supremum  Z  of  processes  of  the  form  (1)  with  Tj^(t)  =  1,  finite  n  and  orthogonally 
invariant  joint  distribution  of  (^i,...,^n)-  In  the  normal  case,  his  results  become 
(Theorem  18.1,  p.  37) 

P{Z  >  a}  ~  ^exp(-a^/2)  (a oo),  (2) 

where  L  =  /o^(I3j=i  and  (Corollary  17.1,  p.  36) 

T  AOO 

>  a}  <  —  exp(-o^/2)  +  /  (27r)"^'^^exp(-x^/2)dx  (c  >  0).  (3) 

27r  Jd 

Johnstone  and  Siegmund  (1989)  consider  processes  of  the  form  (1)  with  Tx{t)  =  1, 
finite  n  and  (Ci»---»^n)  uniformly  distributed  on  the  unit  sphere.  By  making  use 
of  the  well-known  connection  between  the  standaxd  Gaussian  distribution  in  M” 
and  the  uniform  distribution  on  the  unit  sphere  of  K",  we  can  adapt  their  result 
(Theorem  3.3,  p.  190)  to  our  context.  It  turns  out  that  the  resulting  upper  bound 
is  the  expression  (3)  and  then  is  independent  of  n  ! 

Sun  (1993)  investigates  an  asymptotic  expansion  for  the  tail  probabilities  of  the 
maximum  of  smooth  Gaussian  random  fields  with  unit  variance.  In  the  special  case 
of  processes,  her  results  concern  periodic  processes  of  the  form  (1)  with  Tx{t)  =  1. 
For  finite  n.  Sun  obtains  the  asymptotic  formula  (2)  (Theorem  3.1,  p.  40)  as  a 
consequence  of  Weyl’s  Formula  for  the  volume  of  tubes  around  a  manifold  embedded 
in  the  unit  sphere.  For  infinite  n,  (2)  stiU  holds  under  additional  assumptions, 
otherwise  it  becomes  an  upper  bound  (Theorems  3.2  and  3.3,  p.  41). 

Our  approach  is  based  on  the  interpretation  of  the  functions  5j(t),  j  >  1,  as 
a  parameterization  of  a  curve  embedded  in  the  unit  sphere  of  R"  or  of  the  space 
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of  square  summable  sequences.  With  the  canonical  moving  frame  induced  by  this 
parameterization,  we  describe  each  level  manifold  {2:  =  6},  6  e  !R,  of  the  functional 


z  — 


1  ” 

sup  77r  Y.  , 

te/ 


where  {ij  '•  3  >  1}  is  a  realization  of  ;  3  >  1},  cis  an  envelope  of  the  family 
of  hyperplanes  !!Cj=i  =  b  :  t  £  I}.  This  technique  enables  us  to 

express  the  density  fz{b)  of  Z  as  the  canonical  volume  of  {z  =  6},  leading  to  an 
integral  formula  for  fz{b)-  This  method  provides  accurate  lower  and  upper  bounds 
for  fz{b)  and  P{Z  >  c}  (see  Diebolt  and  Posse,  1994,  for  numerical  investigations). 
Moreover,  it  allows  us  to  handle  processes  with  varying  variance,  to  deal  with  the 
case  n  =  00  and  to  manage  the  boundary  effects  with  great  care. 

The  remainder  of  the  paper  is  organized  as  follows.  Our  main  results  are  pre¬ 
sented  in  Section  2  for  n  <  00  and  are  extended  to  the  infinite  case  in  Section  3. 
The  theorems  in  Sections  2  and  3  are  proved  in  Section  4.  Some  related  known 
results  of  differential  geometry  are  briefly  introduced  in  the  Appendix. 

Nota.tion.  Throughout  the  paper,  a.e.  means  almost  every,  refers  to  the  Hilbert 
space  of  square  summable  sequences,  x  =  (xi,X2, ...)  is  an  element  either  of  E" 
or  (x,y)  =  2  <  n  <  00,  t|x||  =  (x,x}^/^  Fect(xi, . .  .,Xd)  and 

Fect-^(xi, . .  .,x<i)  denote  the  linear  subspace  spanned  by  (xi,...,Xii)  and  its  or¬ 
thogonal  respectively,  Gram(xi, . . .,Xd)  is  the  determinant  of  the  matrix  G  with 
entries  =  (a,-,aj)  (note  that  det<J  =  det^(xi,. .  .,Xrf)).  is  the  Gaussian 
measure  on  with  density  V5n(x)  =  (2x)~"/^exp(— ||xlP/2),  ^(x)  =  v?i(x)  and 
=  I-oo  ®y  convention,  <p{oo)  =  0  and  $(cx))  =  1.  C”*(A)  denotes  the 

set  of  functions  A  — ►  E  having  Ath-order  continuous  derivatives  for  k  =  l,...,m. 
The  partial  derivatives  d^^h{x,y)ldx^dy^  where  r  £  is  written  Dkir^x.y). 

The  Jacobian  matrix  of  a  differentiable  mapping  p  :  E"  — E"  is  denoted  Z?p. 


2  Main  results 

Let  X{t)^  t  G  /  =  [0,  r],  be  a  Gaussian  process  with  mean  0  and  variance  <7^(0  >  0, 
of  the  form 

^«)  = 

where  Tx{t)  =  g(t)  =  .  .,gn{t)),  n  >  2,  and  H  =  (6,---,^n)  is  a 

Gaussian  r.v.  with  zero  mean  and  identity  covariance  matrix.  With  the  representa¬ 
tion  (4),  the  covariance  function  rx{ti,t2)  of  X{t)  is  given  by 

Since  crj^t)  =  rxit,t)  ~  ||g(0!P/^x(0i  ilg(01!  =  1  and  g(t)  parameterizes  a  curve 
7  embedded  in  the  unit  sphere  5”~*  in  E’^. 
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In  the  following,  we  assume  that 


(Cl)  the  functions  •rx(<),  5i(t), . .  are  in  C^(I); 

(C2)  Ilg'COII  0  for  all  <  6  J; 

(C3)  {t  :  Cg{t)  =  0}  C  {t  ;  Tx(t)-rjc(t){g'(t),g"(t))/iig'(«)lP  +  r^(f)/llg'(t)iP  >  0} 
where  the  smooth  function  Cg{t)  >  0  defines  the  geodesic  curvature  of  7  at 
the  point  g(f)  (see  Appendix); 

(C4)  whenever  g(t),  g'(t)  and  g{f  )  are  linearly  dependent  for  t'  ^  t,  g{t')  ^  g{t) 
where  g(t)  =  r^(t)g(t)  +  rj^(08'(*)/llg'(0lP- 

Remark.  We  can  relax  condition  (C2)  allowing  a  finite  set  Iq  of  points  such  that 
!lg*(^)!P  =  0  for  t  €  Iq.  (C2)  avoids  additional  technicalities  in  the  proofs  of  our 
results. 

Remark.  If  (Tx(t)  =  is  constant,  (C3)  is  automatically  satisfied,  g(t)  =  Tg{t) 
and  (C4)  means  that  the  curve  7  has  no  self-intersection,  i.e.  t'  ^  t  =>  g{t')  ^  g{t) 
for  all  t,t'  €  (0,T)  if  g(t)  is  T-periodic  and  for  all  t,  t'  6  [0,  T]  otherwise. 

We  are  interested  in  finding  estimates  for  the  distribution  of 

Z  =  sup  A'(t). 
t€/ 

The  key  idea  of  our  approach  is  to  transform  this  problem  into  a  geometric  problem 
concerning  the  standard  Gaussian  measure  of  some  convex  subsets  of  R".  Our  main 
result  is  an  integral  formula  for  the  density  fz  of  Z  which  is  stated  in  Theorem  1. 
The  derivation  of  this  formula  which  is  sketched  below  is  greatly  simplified  if  we 
parameterize  7  with  unit  speed.  This  can  be  done  without  loss  of  generality  when 
condition  (C2)  holds.  Let  us  define  the  Gaussian  process  T(5),  $  €  J,  as 

where  7y(5)  =  tx(A“^(5)),  £(5)  =  g(A~^(5))  and  s  =  A(t)  =  ||g'{u)||  du  defines  a 

unit  speed  parameterization  of  j,  J  —  [0,L]  with  I,  =  I7I  =  A(r).  Then,  we  have 

Z  =  supjf(t)  =  supTi^). 
tel  seJ 

The  covariance  function  of  y(s)  is  given  by  ry(si,S2)  =  rjif(A~^{si),  A~^(s2)). 
Moreover,  ||f(5)||  =  l|f'('S)l|  =  1  for  all  s  e  7.  Note  that  in  terms  of  s  =  A(t), 
(C3)  becomes:  {s  :  Cg{s)  =  0}  C  {5  :  rv(s)  ryl(s)  >  0}. 

For  simplicity,  all  our  results  will  be  expressed  in  terms  of  this  unit  speed  pa¬ 
rameterization.  However,  simple  transformations  give  the  corresponding  formulae 
expressed  in  the  original  parameterization.  In  practice,  only  the  latter  are  used. 
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Our  method  relies  on  the  existence  of  an  orthonormal  moving  frame  (f(s}, 
T(5),Ki(s),  . .  .,K„_2(5))  of  E”  such  that  the  space  tangent  to  5""^  at  f(s)  is 
spanned  by  (T(s),  Ki(s), . .  .jKn-aC-s)).  See  the  Appendix  for  more  details. 

Let  us  consider  the  realization  of  P)  as  (E”,F(E”),/i„)  where  P(E")  de¬ 

notes  the  Borel  cr-field  of  M".  It  follows  that 

y(s,x)  =  — ^(x,f(s)),  XGE", 

Ty(s) 

is  a  realization  of  the  process  K(s)  and 

P{Z  <  a}  =  finiCa),  a  ‘E  E,  (5) 

where 

Ca  =  {x  €  M”  :  supy(s,x)  <  a}. 
s€J 

The  boundaries  dCb  of  Cj,  b  <  a,  partition  Ca-  Indeed,  by  Lemma  6  in  Section  4, 

dCb  =  {x  e  K”  :  supy(s,x)  =  6}. 

This  suggests  that  a  suitable  change  of  variable  wiH  express  tiniCa)  as  an  integral 
over  6  <  a  of  appropriate  superficial  measures  of  dCb- 

Mn(Ca)=  r  MdCb)dh=  r  fz{b)db.  (6) 

oo  •/  — 00 


Such  a  decomposition  can  be  worked  out  basically — for  simplicity,  we  assume 
7  closed  here,  the  case  7  not  closed  can  be  treated  essentially  in  the  same  way — 
because  it  is  possible  (see  Lemma  8)  to  parameterize  dCb  by 

n-2 

P6(s,  u)  =  Ci(h,  s)f(s)  +  c2(6,  s)T(s)  +  Kj(s), 

where  s  £  J,  u  =  (ui,. .  .,Un-2)  €  Db{s),  ci(6,s)  and  C2(b,s)  are  defined  in  terms 
of  b,  Ty(s)  and  its  derivative,  and  Db(s)  is  a  closed  convex  subset  of  We 

show  in  Lemmas  10  and  11  that  the  transformation  p  :  (b,s,u)  pb(s,u)  is  a 
C^-diffeomorphism  from  an  open  subset  of  K”  into  E”.  By  the  change-of-variable 
formula  and  Fubini’s  Theorem,  we  have,  for  all  A  e  B{E"), 


=  /  (pn(p(b,s,u))Gra.m^^^Dp(b,s,u)  dbdsdu 

■/(b,j,u)€p~^(A) 

=  /  /  y’„(p6{s,u))Gram'^^j!)p(6,  s,  u)  dsdu 

JbeM  J(s,u)  €  p^' (AnaCb) 

=  f  tl;b{AndCb)db. 

Jb^m 


Since  Ca  n  dCb  =  dCb  if  b  <  a  and  Ca  n  dCb  =  0  otherwise,  we  obtain  (6). 


Remark.  This  geometric  approach  has  a  probabilistic  interpretation  in  terms  of 
conditional  distributions.  Since  fini-A.)  =  P{H  €  A], 

f^n{A)  =  f  and  P{S  e  A}  =  f  P{S  €  A\  Z  ^  b}fz{b)  db 

JbeM  Vb{oCb)  Jbe* 

for  all  A  €  it  results  that  ipbiA  fl  dCb)/i>b(dCi,)  is  a  regular  version  of 

the  conditional  probability  P{S  ^  A\Z  =  b}.  Similarly,  by  Bayes  formula,  the 
distribution  of  Z  conditioned  on  {E  €  A}  has  density  ipbiA  n  dCb)/ f^n{A). 

Remark.  The  canonical  superficial  measure  on  dCb  induced  by  v’nCx)  is  defined  by 
^^^{AndCb)  ~  I  (pn{Pbis,u))GT&m^^'^Dpbis,u)dsdu 

J{s,^)€p;\AndC^) 

(Berger  and  Gostiaux,  1988,  p.  203).  K  <tk(s)  =  1/r  is  constant,  it  results  from 
Lemma  10  that  -ipbiA  fi  dCb)  =  TwliA  H  dCb)  for  all  A  €  B(R").  Otherwise,  ipb 
appears  as  a  weighted  version  of  with  weight  ry(s): 

‘ipb(AndCb)=  f  ¥3„(P6(s,  «))7y(s)Gram^/^jDpi,(s,  u)ds(iu.  (7) 

7(s,u)€p;'^{Anac4) 

The  above  approach  leads  to  the  following  expression  for  fzib): 

Theorem  1  Under  conditions  (Cl)  —  (C4),  the  density  of  Z  is  given  by 

Mb)  = 

L 

II  TY{s){biTY{s)AT^{s))-uiCgis))ipn(Pib,s,  u))  du^ . .  .dti„-2  ds  +  bzib),  (8) 
0 

where 

.. .  _  r  0  7  closed, 

lTK(0)<,p(6Ty(0))/i„_i(G6(0))  +  Ty(X)9(67y(Z;))^„_i(Gfc(i))  otherwise, 

with 

Db(s)  =  {«  =  (ui,...,Un-2)  sup — )--(p(6,s,u),f(s'))  <  6), 

s'eJ  ^(« ) 

n-2 

p(6,  s,  u)  =  6(Tv(s)f(s)  +  4(s)T(s))  +  ^  «jK_,(s), 

j=i 

Gbil)  =  {u  =  (i;i,...,i;„-i)  e  :  sup  — ^(pt,;(f),f(s'))  <  6}, 

s'^j  ry{s') 

and 

n—2 

PbA^)  =  bTy{l){il)  +  j;n-iT(0  +  5^  t;,K,(/)  /  =  0,  X. 

j=i 
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Remark.  Taylor  expansions  of  order  1  of  (p!i,o(^')*f('*))  ~  bTy{s)  (respectively 
S),.2:,(r;),f(s))  -  bry^s))  around  I,  I  =  0, T,  show  that  Ge,(0)  (respectively  Gi,{L)) 
has  Lebesgue  measure  zero  in  if  7  is  closed. 


Corollary  2  (Upper  bound)  Under  conditions  (Cl)  -  (C4), 
fz{b)  <  M(6)  = 


—  J  Ty(s){TY{s)  +  T^is))exp  |^-y(r^(s)  +  4^(5)) 

+  ^/ Tyis)cg(s)exp  ^-y(4(s)  -|-  4^(5)) j  9  ^ 


$ 


'b{Ty(s)  +  44))' 


Cg{s) 


ds 


f  biryjs)  +  T^{s)y 


ds  +  SM{b) 

(9) 


where 


^M(b)=  I 


0  7  closed, 

7y(0)(^(67V(0))$(64(9))  +  '’y{L)<p{hTY[L)){\  -  $(64(i)))  otherwise. 


Remark.  The  corresponding  expressions  in  terms  of  f  =  A  ^(s)  for  M{b)  and  S\f{b) 
are  obtained  by  replacing  L  by  T,  ds  by  l|g'(t)||(it,  Cg(s)  by  Cg{t)  (see  Appendix), 
ty{s)  by  Tx(f),  r'yis)  by  rlf(()/||g'(i)ll  >>7  -’x(l)(s'(0.S"(<)>/llg’(<)P  + 

Corollary  2  is  a  direct  consequence  of  Lemma  9  in  Section  4  and  can  be  used 
to  derive  a  lower  bound  for  /z(6),  6  >  0.  Indeed,  the  integral  formula  (8)  can  be 
rewritten  as 

fz{b)  =-^J n'(5)(Ty(s)  +  4(s))exp  +  4^(«)))  (^n-2{Db{s))  ds 

n  ' 


L 

1  /■  /  X 


Cg{s)exp  \ 


b\ 


(t^(s)  +  4^(^)) )  /  “iSPn-2(w)  dui...  dun-2  ds+6z{b). 


(10) 

By  definition  of  Db{s)  and  relation  (5),  y.n-2{Db{s))  and  Uiy)„_2(M)dt£  can  be 
interpreted  in  terms  of  /’(sup,<gj  W^»(50  <  b),  where  1^5(5')  is  a  suitable  Gaussian 
process  of  the  form  (4)  for  a.e.  s  £  J.  Therefore,  we  can  make  use  of  Theorem  1  and 
Corollary  2  to  provide  upper  bounds  for  P{supy^j  W^s(s')  >  a)  and  the  absolute 
value  of  the  second  term  in  the  right-hand  side  of  (10).  This  approach  requires  that 
the  process  ITs(s')  satisfies  conditions  (C1)-(C4).  The  resulting  assumptions  are 
still  mild  but  very  technical.  For  simplicity,  we  give  sufficient  conditions  for  n  >  5: 


(C5)  the  function  Tx{t)  is  in  C^(/)  and  the  functions  gi(t), . .  -ignit)  are  in  C%I); 
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(C6)  for  aU  i  €  I,  Tx{t)  -  r^0)(g'(0,g"(0}/llg'(0lP  +  r’i(t)n\g'{t)f  #  0; 

{C7)  for  all  f  €  I,  =  rx{t')  -  rxCiXsCO.gCO)  -  (0(g'(<).g(«'))/l|g'(<)li^ 

>  0  for  all  t'  t'  G  /; 

(CS)  for  a.e.  t  €  /,  the  vectors  g(f),  g'(t),  g(t'),  g'(t')  and  g"(t')  are  linearly 
independent  for  all  t'  ^  t,  t'  €  /; 

(C9)  for  a.e.  t  €  I,  the  vectors  g(t),  g'(t),  g(t'),  g’{t')  and  g(t")  are  linearly 
independent  for  aU  t'  ^  t,  t"  t,  t"  /  t'  and  t\t*'  £  I. 

Remark.  In  terms  of  s  =  X{t),  (C6)  and  (C7)  say  that  for  all  s  G  J,  7y(s)+ry  (s)  ^ 
0  and  /3s(s')  =  Ty{s')  -  Ty(s)(f(s),f(5'))  —  ry(s)(T(5),  f(s'))  >  0  for  all  s'  5^  s. 
s'  G  J-  (C7)  implies  that  Ty(s)  +  ry(s)  >  0  for  aU  s  G  ^  since  3s{s  +  k)  = 
(tv(s)  +  TY(s))h^f2  +  o(h^)  as  ft  — >  0. 

Remark.  Under  (C8),  l|g'(f)|!  ^  0  and  c^(t)  >  0  for  all  t  G  /,  and  (C2)  and  (C3) 
applied  to  the  process  Wj(s')  hold  for  a.e.  s  €  J-  Under  (C9),  (C4)  applied  to  X{t) 
and  1^5(5')  for  a.e.  s  G  J  holds.  (C7)  implies  (C4).  If  axit)  =  l/r  is  constant,  the 
reverse  implication  is  true. 

More  precisely,  let  us  define  the  functions 


o,(/)  =  (l-{f(s),f(/))=-(T(s),f(»'))")''*, 

(11) 

7-5(5')  =  3,($')/as(s'), 

(12) 

T}{s)  =  inf  T,(s'), 

(13) 

r5(5'i,s'2)  =  ar^(-Si)a7H4) 

X  ((f(si),f(s'2))  -  (f(s),f(s;)){f(s),f(s'))  -  (T{s),f(si))(T(s),f(si)))  (14) 
if  s'  s  and  s^  5^  s,  r,(5,s'2)  =  (f(s)  +  f"(s),f(s^))/(c5(s)a5(s'2))  and  r,(s,s)  =  1. 

By  condition  (C9),  a,(s')  =  (i:;)=i^(K,(s),f(s'))2)V2  >  0  for  all  s'  5^  s,  s,  s'  G  J. 
By  condition  (C5),  r,(si,S2)  has  continuous  partial  derivatives  Z?w^j(>s'i,S2)  for 
0  <  k,l  <  4.  Moreover,  r,(s')  >  0  and  rj{s)  >  0  for  all  5,  s'  G  J. 

For  a.e.  s  G  J,  W',(s')  =  Tj“^(s')(0,ks(s')),  where  fi  =  (u)i, . .  .,u;„_2)  is 
a  Gaussian  r.v.  with  mean  zero  and  identity  covariance  matrix  and  ks(s')  = 
(^s.iC^'),  - .  fcj,„_2(s'))  parameterizes  a  curve  7,  on  the  unit  sphere  of  More¬ 

over,  r,(si,S2)  =  (k,(si),ks(52)}.  This  curve  is  the  normalized  orthogonal  projec¬ 
tion  of  7  on  Fec<-‘-(f(s),T(s))  and  consequently  s'  -+  ks(s')  is  not  unit  speed.  In 
order  to  apply  Corollary  2  to  Z,  =  sup,fgj  W',(s'),  it  remains  to  determine  l|k'(s')|j 
and  the  geodesic  curvature  €3,5(5')  of  7*. 

Let  us  denote /?«, 5(5')  =  Dkir,{s'^,s'2)\s[=:>'^=s'-  We  have  !lk;(5')|i  =  p\y,{s')  and 

=  {pnA^')p22As')  -  Pi2,,(s'))//>ii,,(«')  -  1- 

Theorem  3  (Lower  bound)  Under  conditions  (C5)  -  (C9)  and  for  n  >  5, 
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fz{b)  >  Tn{b)  = 


^  J  7y(s)(Tv(5)  +  r{J(s))exp^-y(r^(s)  +  rf2(^s))^  ^  -  J  M,{b')db^  ds 

J 7V(«)c3(s)exp  j(r^(5)  + 4^(s))^  77-*(5)A/j(6)(i5.  (15) 


_  _L 

27r 

for  all  b  >  0,  where 
Ms 


,'(»)  =  |^/r.(»')C,(»')«P  {-J  (r;(»')+  ♦  (^)  pYZis'W.-' 

+  ^/r.(^0P,/*0exp(-^(r;(/)+£i^))9(^)p!C,(»0*'+<M.(i), 

{0  7  closed, 

r,(0)¥>(6r,(0))$(j|^)  +r,(i)¥>{Mi))  (l  "  otherwise, 

and 

(sis')  =  rs(s')  -  ri(s')pi2,5(5')/Pnfi(«')  +  ■^"(5')/pn.*(«')- 


Remark.  The  corresponding  expressions  in  terms  of  t  =  A~^(s)  are  obtained 
by  using  the  same  transformations  as  in  Corollary  2  and  by  replacing  r,(s')  by 
n^t')  =  A(l')/«((<')  "kere  offC)  =  1  -  (g(().g(''))“  -  (g'(<).g(<'))^/llg'(l)ll'-  ’■;(«') 
by  »-,'((')/llg'(*')ll,  '■."(»')  by  -T;((')(g'(('),g”((')>/i|g'(<')|P  +  -•;'((')/llg'(‘')ll"  and 
r,(si,S2)  by  ^<(^1,^2)  with  similar  transformation  rules  for  Pki.a{s'). 

Together  with  Corollary  2,  Theorem  3  enables  us  to  investigate  fz{b)  and 
P{Z  >  a}  with  a  new  accuracy.  Both  upper  and  lower  bounds  for  the  tail  dis¬ 
tribution  of  the  maximum  of  large  classes  of  Gaussian  processes  have  been  obtained 
by  Samorodnitsky  (1991)  and  Berman  and  Kono  (1989)  among  others.  However, 
these  results  are  asymptotic  and  involve  unknown  constants.  The  next  theorem 
states  the  efficacy  of  our  bounds.  It  shows  that,  in  general,  they  are  asymptotic 
to  fzib)  as  b  —>■  00.  In  addition,  simple  asymptotic  formulae  for  fz(b)  are  made 
explicit. 


Theorem  4  Under  the  conditions  of  Corollary  2  for  M{b)  and  of  Theorem  3  for 
m{b),  and  consequently  for  fzib), 

fzib)  =  e(6)(l-bi2(6)) 

M(b)  <  eib)il  +  RMib)),  i?Af(h)>0, 
mib)  >  eib)il- Rmib)),  -R^h)  >  0, 

as  6  — ♦  00,  where 
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(t)  i/7  is  closed,  or  7  is  not  closed  and  min(Ty(0).  Ty(I))  >  ry(s). 

e(6)  =  ^  ^  rY{s){Ty{s)  +  TY{s))exp  +  rf  2(5))^  ds  (16) 

and  R{b)  =  Rm^^)  =  Rm{b)  -  O{ip{b0))  for  some  9  >  0; 

{ii}  if ')  is  not  closed  and  Ty{s)  =  1  is  constant, 

e{b)  =  ~Lex'p{-9l2), 

R(b)  =  RM(b)  =  and  RM{b)  =  0{!fi(b6))  for  some  9  >  0,' 

(Hi)  if  is  not  closed,  Ty(0)  =  r  <  7y(s)  for  a//  s  ^  0  and  q  =  mf{fc  >  1  : 
ry^^(O)  =  >  0}  w  assumed  finite,  then 

1.  e(b)  =  o^“‘/’exp(-6*r^/2)(27r)“^r‘~^'''’(r(^'5)“^/«g“^r(g“^)r(g+l)*^^  and 
R(b)  =  RM{b)  =  Rm(b)  =  0(1)  ifq>  3; 

e(6)  =  exp(-6VV2)(2;r)-i/2r(2-Hr")-'/'(T+r")'/'*+l),  Rib)  =  Rmib)  = 
0(1)  and  Rrr^ib)  =  1  -  (2-H^")~'/'(^  +  +  1)''  -  o(l)  if  q  =  2; 

5.  e(6)  =  Suib),  Rib)  =  o(l),  iZ^Cfr)  =  Oi<pib9))  and  il^(6)  =  1  -  0(<^(6^)) 
/or  some  9  >  0  if  q  =  1. 

Remark.  The  remainders  Rib),  R\fib)  and  Rmib)  are  of  order  o(l)  whenever  we 
have  used  Laplace  approximations  (De  Bruijn,  1962,  p.  65)  in  (16)  to  determine  the 
asymptotic  formula  for  fzib)  (see  proof  of  Theorem  4).  Such  occurrences  do  not 
reflect  an  inaccuracy  of  our  bounds. 

Theorem  4  shows  that,  in  general,  fzib)  can  be  approximated  with  (unexpected) 
high  accuracy  by  the  simple  expressions  e(6).  It  also  reveals  that  when  the  process 
XU)  (T(s))  reaches  its  maximum  at  the  boundaries  with  high  probability,  the 
main  contribution  to  the  density  fzib)  is  given  by  the  additional  term  Szib).  This 
phenomenon  affects  the  good  behavior  of  m(b)  since  we  have  chosen  to  take  0  as 
lower  bound  for  Szib)  for  sake  of  brevity.  However,  it  is  possible  to  improve  m(6)  by 
introducing  a  term  0  <  Smib)  <  Szib)  which  corrects  this  imperfection  (see  Diebolt 
and  Posse,  1994). 


3  Extension  to  series  representation 

The  results  of  Section  2  can  be  extended  to  the  maximum  Z  of  Gaussian  processes 
of  the  form 

Xii)  =  rJ^^it)Uit)  t€/  =  [0,T], 
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where  Tx{t)  >  0,  Uit)  is  a  centered  Gaussian  process  with  anil  variance  and  finite 
covariance  function 

OC 

(g(<i),g(f2)}  =  Yl3AU)9:if2)- 

j=i 

The  functions  g(/)  =  (si(Gi?2(0!  •  •  •)  parameterize  a  curve  'y  embedded  in  the  unit 
sphere  of  £^.  Let  us  assume  that 

(Dl)  the  function  rx(t)  is  in  C^{I)  and  rL:(ti,t2)  bas  continuous  partial  derivatives 
DkiTvihAi)  for  0  <  A:,  /  <  4. 

Let  us  denote  Pki,A^)  =  Dkiru{h,t2)\ti-t2-f  We  also  assume  that 
(D2)  Pii.uCO  7^  0  for  all  t  e  /; 

(D3)  Tx{t)  -  Tx(0pi2.«(0/Pn?«(0  fi  ^x(0/pii.u(0  >  0  for  all  t  €  /; 

(D4)  the  curve  7  parameterized  by  g(t)  =  rx(t)g(t)  +  Tx(0g^(0/Pn.u(0  bas  no 
self-intersection. 

Under  these  conditions,  we  show  the  existence  of  a  density  for  Z  and  give  an  upper 
bound.  With  the  following  assumptions  we  also  derive  a  lower  bound.  Note  that 
condition  (D5)  can  be  weakened  since  it  is  imposed  by  the  sufficient  conditions  (D7) 
and  (D8)  (see  Section  2). 

(D5)  The  function  Tx{t)  is  in  C*{I)  and  ruih,  12)  has  continuous  partial  derivatives 
DkiTvih-ih)  for  0  <  A;,*  <  14; 

(D6)  for  all  t  €  /,  =  Tx{t')  -  Tx{t)Tv{tA')  -  T‘x{i)DiQru{t,t')l PuM)  >  0 

for  all  t'  ^  t,  t'  €  /; 

(D7)  for  a.e.  t  e  I,  the  sequences  g(t),  g'(t),  g(t')>  ^.nd  g"(t')  are  linearly 

independent  for  all  t'  7^  t,  t'  €  I; 

(D8)  for  a.e.  t  £  I,  the  sequences  g(t),  g'(t),  g(t')i  s'{f')  a^nd  g(<")  are  linearly 
•ndependent  for  all  t'  ^  t,  x"  ^  t,  t"  ^  t'  and  t',  t"  £  I. 

The  length  and  geodesic  curvature  of  7  are  given  by  L  -  Jq  and  cj{t)  - 

(Pn,«(0p22,«(0  -  Pl2,u(0)/Pll.u(<)  -  1- 

As  in  Section  2,  7  has  a  unit  speed  parameterization  s  =  X{t)  =  Jo  Pu\it')dt' 
under  (D2).  Therefore,  Z  =  sup,g/A’{t)  =  sup,g jy(s)  where  Y{s)  =  r(7^(s)V'(s) 
with  7y(5)  =  rx(A“^(s)),  F(s)  =  Z7(A“’(s))  a  Gaussian  process  of  variance  one 
and  covariance  function  rv'(si,S2)  =  (f(si),f(s2))  =  and  f(s)  = 

g('’'"H'S))-  We  have  c^Cs)  =  /?22,v(s)  -  1  where  PkiAs)  =  T>i;rv-(si,S2)U,=:,j=,. 
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Finally,  the  functions  as(s').  >7(s),  r,{s{,S2),  PkiA^')^  Cp.,(s')  and  Cs(s')  de¬ 

fined  in  Section  2  are  also  well-defined  in  the  present  context  and  we  can  show  the 
following  result. 

Theorem  5  (i)  C'^nder  conditions  (Dl)  -  (D4),  Z  has  a  density  fz{b)  and  fzib)  < 
M{b)  for  all  b,  with  M{b)  given  by  (9); 

(ii)  Under  conditions  (D3),(D5)  —  (D8),  fzib)  >  m(b)  for  all  b  >  U,  with  m(b) 
given  (15). 

In  addition,  Theorem  4  still  holds. 


4  Proofs  of  results  of  Sections  2  and  3 

4.1  Proof  of  Theorem  1 

We  give  a  detailed  proof  for  the  case  7  closed  and  sketch  the  straightforward  adap¬ 
tation  for  the  other  case. 

We  need  a  workable  description  of  the  boundary  dCa  of  Ca,  for  a  €  R.  Note 
first  that  Ca  is  the  intersection  of  the  closed  half  spaces  {x  €  R”  :  y(s,x)  <  a}, 
s  £  J.  These  half  spaces  have  hyperplane  boundaries  ^ven  by 

Fr(s;  c)  =  {x  €  R"  :  F(s,x)  =  a},  $  £  J.  (17) 

Lemma  6  shows  that  the  surface  dCa  is  closely  related  to  the  hypersurface  enveloped 
by  the  hyperplanes  H{s-,a),  s  £  J. 

Lemma  6  IfCa^  9,  dCa  =  {x  €  R"  :  sup,gy  y(s,x)  =  a}. 

Proof.  dCa  D  {x  €  R"  ;  sup,gjy(s,x)  =  a}.  Let  x  G  R"  satisfy  sup,gjy(s,x)  = 
a.  Then  x  £  Ca.  Let  us  prove  that  it  does  not  intersect  int{Ca).  Since  the  function 
s  €  J  — ♦  y(s,  x)  is  continuous  on  the  compact  interval  J,  it  reaches  its  maximum. 
Therefore,  there  exists  sq  £  J  such  that  sup,£jy(s,x)  =  y(so,x)  =  a.  Let  c  >  0 
be  an  arbitrary  small  positive  real  number.  In  the  open  ball  i5(x,2e)  with  center  x 
and  radius  2c,  we  can  pick  y  =  x  -I-  cf(so).  Then  y(so,y)  =  a  +  7^;^^  >  a-  Hence 
y  ^  Ca  and  x  £  dCa- 

dCa  C  {x  €  :  sup,gjy(s,x)  =  a}.  We  will  show  that  if  sup,gjy(s,x)  <  a, 

then  X  €  int(Ca).  Since  the  function  y(s,x)  is  continuous,  there  exists  sq  £  J  such 
that  supjgj  1  j,x)  =  y(so5x)  =  b  <  a.  Then,  for  any  sufficiently  small  c  >  0  and 
for  each  y  £  i?(x,  c),  we  have 

supy(s,y)  <  6-f  sup-4— (y-x,f(s)) 
aeJ  Ty(«) 
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< 


1 


< 


b  +  sup 


eJ  Ms) 


sup  sup  |{z,f(5))| 
2eS(0,£) 


b  +  sup 


1 


sej  Ty{s) 


(  <  a. 


0 


Lemma  7  Ca  =  Sa  flCa.  where  Ea  =  U^gjIInCs),  with  Ilats)  denoting  the  affine 
subspace  of  dimension  n  —  2  ofW'  defined  by  the  equations 

r  (x,f(s))  =  <i7y(s) 

1(x,T(5))  =  aMs). 


Proof.  dCa  C  Eo  n  Co-  From  Lemma  6,  it  follows  that  for  each  x  6  dCa,  there 
exiSts  So  €  J  such  that  sup,gj  y(s,x)  =  y(so,x)  =  a,  i.e. 

{x,f(so))  =  a7v(so)-  (19) 


At  such  a  point  so,  the  function  s  — ► 

dYis,x)  _  r^(so) 
ds  T-^C-so) 


y (s,  x)  reaches  its  maximum,  implying  that 
(x,f(so))  +  ;j^(x,T(so))  =  0.  (20) 


Plugging  (19)  into  (20),  we  obtain  (18). 

dCa  3  Sa  n  Cfl.  Let  X  €  Ea  HCq.  Then,  sup,gjy(s,x)  <  a  (since  x  €  Co) 
and  there  exists  sq  6  J  such  that  y(so,x)  =  a  (since  x  €  Eo).  This  implies  that 
sup3gjy(s,x)  =  a,  which  in  turn  implies  that  x  €  dCa,  by  Lemma  6.  0 


Lemma  8  (i)  The  hypersurface  Ea  con  be  parameterized  by 

n-2 

Pa(s,  u)  =  a(Ty(s)f(s)  +  4(s)T(s))  +  ^  tijKj(s),  (21) 

J=1 

with  s  £  J,  u  =  (ui , . . . ,  Un-2)  € 

(n)  The  hypersurface  dC^  can  be  parameterized  by  pa(5,  u),  with  s  £  J  and  u  € 
Z?a('S)?  where  Da{s)  is  the  closed  convex  subset  o/K"~^  (possibly  empty)  defined  by 
the  set  of  inequalities  sup,»gy  Y (s',  Pa(s, «))  <  a. 


Proof,  (i)  Since  Ea  =  U,6jna(s),  it  suffices  to  prove  that  for  each  given  so  €  J, 
Pa(so,M),  u  £  parameterizes  the  affine  subspace  no(so).  First,  Pa(so,«)  € 

na(so),  since  it  satisfies  (18).  Second,  the  function  u  -*  Pa(so,  w)  is  an  affine  map 
from  to  Ilalso).  Since  the  dimension  of  the  affine  subspace  Ilalso)  of  1”  is 
n-2,  it  remains  to  prove  that  u  — ►  pa(5o,w)  is  an  affine  bijection.  This  is  true  since 
Ki(so),  ■ . . ,  K„_2(so)  are  linearly  independent. 

(it)  This  foOows  directly  from  (i)  and  Lemma  7.  0 
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Lemma  9  Let  sq  £  J  be  given. 

(i)  Ifu€  Da(so)  7^  0.  then  d(a,so)  -  uiCg{so)  >  0, 

(ii)  if  u  £  int{Dc:iso))  ^  0  and  Cg{so)  >  0,  then  d{a,SQ)  -  UiCg{so)  >  0,  where 
d{a,s)  -  a{TY{s)  + 

Proof,  (i)  For  each  fixed  u  £  and  3ot  the  function  hu.soCs)  =  F  (s.ps{5C)  «))!•?€ 

J,  is  twice  differentiable  and  =  0.  Furthermore,  since  f"(s)  =  Cj{s)K(s)  ~ 

f(s),  for  all  s  €  ~  -(d(a,so}  -  uiC5(5o))/7y(5o)-  If  u  €  -Da(-So)  #  0. 

/iu,so(^)  reaches  its  maximum  value  at  s  =  sq,  implying  that  h'^  ,^(so)  <  0. 

(ii)  Suppose  that  u  £  int(Daiso))  ^  0.  Let  us  show  by  contradiction  that 

K, 3oi^o)  <  0.  Otherwise,  we  would  have  /»u,,o(so)  =  0  by  (i).  If  h"g^{so)  =  0,  since 

Cj(5o)  >  0  and  u  £  int(Da(sQ)),  we  can  pick  i;  £  Da{so)  (close  enough  to  u)  such 
that  /i",jo(so)  >  0  (by  taking  vi  >  «i),  which  contradicts  (i).  0 

Remark.  If  0^(50)  =  0,  the  inequality  h*^  j^{so)  <  0  can  hold  iff  d(a,so)  >  0. 
This  shows  that  in  this  case  L)a(so)  =  0  for  d(a,So)  <  0.  If  d(a,so)  >  0,  then 
K3o(^o)  =  -d(a,so)/Ty(so)  <  0  for  all  u  € 

Let  us  define  the  function 

p(6,s,u)  =  pj,(s,u),  b€M,  s£j,  tz  =  (ui,,.., tin-2)  e 

This  function  maps  Va  =  {(6,s,u)  :  b  <  a,  s  £  J  and  u  £  Di{s)}  C  HR”  onto 
Ca  C  K”  since  Ui,<o5Ci  =  Ca-  Moreover,  the  restriction  of  p(f»,s,«)  to  the  open 
subset  =  {(f>,5,ti)  \  b  <  a,b  ^  Q,s  £  int{J)  and  u  £  int{Dbis))}  C  K  maps 
onto  a  subset  C®  of  Co.  In  the  following,  we  will  assume  that  a  is  such  that  ^  0 
(by  convention,  if  tnt(l?6(s))  =  0,  then  (b,s,u)  ^  F®). 

Lemma  10  (t)  |  det  £>p(6,  s,u)|  =  Tv(5)ld(6,s)  -  i/iCj(5)|. 

(ii)  The  function  p(h,  s,ii)  is  a  local  -diffeomorphism  from  onto  C°  and  C® 
is  an  open  subset  o/R”. 

Proof,  (i)  It  suffices  to  prove  that  Gram  Dp  =  ry(s)(<Z(6,s)  —  tiiCj(s))^.  This 
is  true  since  dp/db  =  Ty(s)f(5)  +  ry(s)T(s),  dp/ds  =  (d(b,s)  -  tiiCj(s))T(s)  + 
ej(s)Kj(s),  dpiduj  =  Kj,  j  =  -  2  and  the  functions  e_j(s),  j  = 

1, . . . ,  n  —  2,  need  not  be  known  since  they  disappear  in  the  evaluation  of  Gram  Dp. 

(ii)  It  follows  from  (z)  and  Lemma  9  that  for  (b,  s,  u)  £  F®,  det  Dp(b,  s,  u)  /  0. 
Hence,  the  restriction  of  p  to  F°  is  a  local  C^-diffeomorphism  from  F®  onto  it.s 
image  C®  and  C®  is  an  open  subset  of  M”. 

Lemma  11  p  is  a  one-to-one  mapping  from  F°  onto  C°. 

Proof.  Since  the  hypersurfaces  dCb,  b  <  a,  partition  Ca,  it  suffices  to  examine  the 
case  where 

p(6,so,w)  =  P(6,5i,v),  (22) 
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with  {b,so,u)  and  {b,si,v)  in  Moreover,  if  sq  ~  Si,  u  =  c  since  Kj(£o). 
j  =  1, . . . ,  n  -  2  are  linearly  independent. 

It  follows  from  (22)  with  sq  .si,  soj'Si  6  int(J)  that  }''{si,p(b,SQ,u))  =  b.  If 
there  exists  y  (1  <  j  <  n-2)  such  that  (Kj{so)5  a  perturbation  argument 

similar  to  the  one  used  in  the  proof  of  Lemma  9  {ii)  shows  that  Y  (si,p(b,  sq,  u))  <  b. 
since  u  e  int(Di,(so)).  If  (Kj(so)5f(-si))  =  (K_,(si),f(so))  =  0  for  aU  j  (1  <  j  < 
n  -  2),  then  f(si)  e  V"ect(f(so),T(so)),  f(so)  €  Vect(f(si),T(si))  and 

6(7v(so)^(5o)  +  ^y(^o)T(so))  = 

which  is  impossible  since  b  ^  0  and  condition  (C4)  holds.  0 

Lemma  12 

Mn(C'a)  =  /  ry(s)(d(b,  s)  -  uiCg(s))<pn(p(b,  s,  u))  dui...  durt-2  ds  db.  (23) 
Jvo 


Proof.  According  to  Lemmas  10  and  11,  the  function  p  is  a  C^-diffeomorphism  from 
V°  onto  C°.  Moreover,  according  to  Lemmas  9  and  10,  detDp(b,  s.u)  >  0  for  all 
(b,s,u)  €  V'®.  Then,  (23)  results  from  the  change-of-variable  x  =  p(6,  s,  u)  applied 
to  the  integral  /iin(C°)  =  f^o  ‘Pn(x}dx.  C> 

Since  Pi(s)  is  a  convex  of  the  Lebesgue  measure  in  of  dDb(s)  is  zero. 
Therefore,  the  Lebesgue  measure  of  Va\^<?  is  zero  and  we  can  replace  by  Va  in 
(23).  Moreover,  Ca\Ca  has  Lebesgue  measure  zero  since  Ca\C°  =  p(i4)\p(L’°)  C 
p(V’o\L),°)  and  p(6,5,u)  is  a  C^-function  from  M"  to  R”.  Consequently,  /x„(Ca)  = 
Pn(C'o)  which,  with  (5)  and  (6),  concludes  the  proof  of  Theorem  1.  0 

Suppose  now  that  7  is  not  closed.  Lemma  6  still  holds.  In  Lemma  7,  we  have 
to  replace  Ea  by  E^  =  (Ujginf(J)na(s))U  H(0,a)U  H{L,a)  where  n{La)  is  defined 
in  (17).  Then  dCa  can  be  partitioned  as  dCa  =  dCaMt  U  dCa.o  U  dCa,L  where 
dCa,int  =  (Usg,„t(j)no('S))  n  Ca,  dCa,i  =  H{l,a)n  Co,  1  =  0,i.  Lemmas  8-12  can 
be  applied  without  modification  to  Ca,int  ~  ^b<adCa,int-  The  additional  term  6z{b) 
is  obtained  from  the  boundaries  Co,/  =  Ub<a9Ca,h  1  —  0,1.  Indeed,  for  I  =  Q.L. 
dCa,i  can  be  parameterized  by  P6,/('y)  =  6Ty(/)f(/)+  Vn-iT(/)  +  VjKj(l)  with 
V  €  Gb{l)  =  {v  €  :  supj,gj7y(s')~^(p6,/(u),f(s'))  <  f>}.  Since  Gram  Dpb.i  = 

4(0,  we  have  il:b{dCb,i)  =  '^n{p6,/(v))Ty(/)du.  0 

4.2  Proof  of  Theorem  3 

Let  J*  be  a  subset  of  J  such  that  J\J*  has  Lebesgue  measure  zero  and  for  each 
s  €  J*,  f(5),  T(s),  f(5'),  T(s')  and  N(s')  are  linearly  independent  for  all  s'  ^  s, 
s'  G  J,  and  £(«),  T(s),  f(s'),  T(s')  and  f(5")  are  linearly  independent  for  all  s'  s, 
s"  ^  s,  s"  5*^  s'  and  s',  s"  G  J. 
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For  each  s  €  </“,  ^in-2(Db{s))  in  formula  (10)  can  be  interpreted  as  P{Zs  <  b) 
where  Z,  =  sup^/gj  W',(s'), 


“i'»(-s)/(n'(-s) +  >■?(«)) 


s'  ^  ^ 
s'  =  s, 


with  n  =  (uj, . .  .,a;„_2)  a  Gaussian  r.v.  with  zero  mean  and  identity  covariance 
matrix.  With  this  definition,  Ws{s')  is  defined  and  continuous  on  J,  VarlW^is')}  = 
al{s')/0^{s')  for  s'  ^  s  and  Var(Ws(s))  =  Cj(s)/(7y(s)  +  ryJ(s))2.  Moreover,  Wj(s') 
is  of  the  form  (4)  and  satisfies  conditions  (Cl)-(C4):  1^,(5')  =  k3(s'))- 

where  Ta{s')  =  {Var(Ws{s')))-^^'^  and  k,(s')  =  iks,ii$'), A:i,„_2(s'))  with  ks,j{s') 
=  (Kj(s),f(s'))/as(s')  for  s'  «,  .ta,i(s)  1  and  k^^j{s)  =  0,  i  >  2.  Therefore,  we 
can  apply  Corollary  2  to  Zs  to  obtain  an  upper  bound  for  1  -  ^n_2(Z)(,(s)). 

The  above  interpretation  can  also  be  used  to  derive  an  upper  bound  for  the 
second  term  in  formula  (10).  Indeed,  by  Stokes’  Theorem  (Berger  and  Gostiaux, 
1988,  page  195) 


J  Ui<Pn-2iu)  dUi...  dUn-2 

— 

^  (Pn-2(u)  dU2  A  ...  A  dUn-2 

Db{s) 

SDiia) 

<  J  'Pn-2  dv, 

dD^(s) 


where  dV  denotes  the  canonical  volume  element  of  the  manifold  5Di,(s)  (Berger  and 
Gostiaux,  1988,  page  203).  A  straightforward  adaptation  of  Lemmas  6-9  yields 

/  <fin-2dV  =  f  SC'n-2(P6,»(5',i^))Gram^''^I>P6,,(s',u)d5'dv, 


where  s'  €  J,  u  €  R”  ^  and  (s',  v)  -+  P6,3(s',  u)  defines  a  parameterization  of  dDbis) 
analogous  to  (21).  From  (6)  and  (7)  applied  to  the  process  Ws(s'),  it  follows  that 


L 


dDi{,) 


V’n-Z  dV  < 


inf  tJs') 


rJz.W  =  n-\s)fz,{b), 


where  fz,{b)  density  of  Z,. 


0 


4.3  Proof  of  Theorem  4 


(t)  Assume  7  closed  and  let  denote 
L 


eiW  =  ^  /  w(5)(Ty(s)  +  r^(s))exp  ^-y(r^(s)  +  4^(5))^ 


By  (9),M(6)  =  ei(5)  + 
L 


J  ry{s)iTY{s)  +  T^{s))exp  f-^(4(5)  +  4^(5))  j  H  ds. 


(®) 
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where  0  <  H{x)  =  x~^ip{x)  —  (1  —  $)(ar)  <  x~^'^(x)  for  all  x  >  0  and  H{oc)  =  0. 
Since  Cg(s)  is  continuous  and  7y(s)  +  ry(5)  is  continuous  and  positive  on  J .  it  follows 
that  9i  =  infsgj{Ty(s)  +  •ry(s))/cg(s)  >  0  and 

M(6)-e,(6)<Ci(6)(60i)-V(i^i). 


By  (15),  m(6)  =  e\{b)- A- B,  where  A  involves  Ms{b')dh'  and  B  involves  Mj(6). 
Therefore,  it  suffices  to  examine  the  asymptotic  behavior  of  Since  t){s)  >  0. 


Ms{b)  <  ^exp  J  r,(s'){bCsis')  +  c,,s{s’)(2n)-^^^)pl{lis')ds'. 


Using  Taylor  expansions  of  sufficient  order,  it  can  be  shown  that  the  functions 
{$,$')  ^  J  y.  J  ^  Ts{s'),t'Xs'),t"{s*),  Pki,s{s'),  0  <  fc,/  <  2,  and  Cg^^is')  are  con¬ 
tinuous.  Therefore,  their  supremum  over  J  x  J  is  finite.  Moreover,  by  the  conti¬ 
nuity  and  positivity  of  rj(s')  as  a  function  of  (s,s')  6  J  x  J,  ^2  =  infagj  r}{s)  = 
inf(i,4')g7xJ'rs('SO  >  0-  follows  that  <  Cibip(b62)  for  all  s  €  J  and  all 

6  >  0,  hence  Ms(b')db'  <  C2‘p(b02).  Consequently,  A  <  61(6)039(6^2)  and 
B  <  61(6)049(6^2)-  From  the  asymptotic  behavior  of  M{b)  and  m(6),  it  results 
that  ei(6)  =  e(6). 

If  7  is  not  closed  and  min(Ty(0),Ty(T))  >  inf^gj  Ty(5), 


9(6ry(0)) 
61  (6) 


Cs  hy(s){TYis)  +  T^{s))  exp 


(r^(0)  -  4{s)  -  42(5)) 


For  €  >  0  small  enough,  the  subset  =  {s  £  J  :  ry(0)  —  ry(s)  —  ry^(s)  >  6^} 
of  J  contains  a  nonempty  interval  (around  a  global  minimum  of  ry(s)  -f  ry^(s)). 
Therefore,  it  has  positive  Lebesgue  measure.  For  such  an  e  >  0, 

Ty(s)(ry(s)  +  4(s))exp  ^y(4(0)  -  4(s)  -  4^(5))  j  ds 

>  exp(6V/2)  /  Ty(s)(Ty(s)  -1-  T^{s))ds  =  Cev^~\b€). 

J  Jf 

Hence,  9(6ry(0))  =  61(6)0(9(66*3)),  i.e.  6a/(^>)  =  61(6)0(9(6^3)).  To  complete  the 
proof,  it  remains  to  show  that  sup,£j6a/,(6)  =  61(6)0(9(6^4)).  This  follows  from 
the  continuity  and  positivity  of  s  — ^s(O)  and  r,(X). 

(if)  Straightforward. 

(iii)  First,  note  that  ^(0)  >  0.  Since  rn(6)-  ei(6)  =  61(6)0(9(6^5))  by  (i), 

it  suffices  to  study  the  asymptotic  behavior  of  Ci(6)  and  6m(&)-  The  result  follows 
directly  from  the  Laplaiie  approximation  (De  Bruijn,  1962,  p.  65)  of  ei(6)  using  a 
Taylor  expansion  of  order  q  of  4(s)  -f  4^('®)  ^-round  s  =  0.  <> 


/ 
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4.4  Proof  of  Theorem  5 


Lemma  13  shows  that  the  process  X{t)  admits  the  representation  (1).  Therefore,  by 
truncation  and  renormalization,  we  can  construct  a  sequence  of  Gaussian  processes 
X„(i)  of  the  form  (4)  which  converge  to  X(t).  Moreover,  under  (D1)-(D4),  Xn(t) 
satisfies  (Cl)-(C4)  for  all  n  sufficiently  large.  Similarly,  under  (D3),  (D5)-(D8), 
Xn(t)  satisfies  (C5)-{C9)  for  n  large  enough.  Hence,  the  density  fzni^)  of  = 
suPig/ Jtn(f)  has  an  upper  bound  M„(6)  of  the  form  (9)  and  a  lower  bound  mn{b) 
of  the  form  (15). 

We  will  show  that  Mn(fe)  — >  M{b)  for  all  b,  m„(6)  m{h)  for  all  6  >  0,  and 
the  sequence  {fzn}  is  weakly  relatively  compact  in  With  an  equality  due  to 

Dmitrovskii  (Lifshits,  1986),  this  implies  that  Z  has  a  density  fz{b)  which  is  the 
limit  of  fz„{b)  and  satisfies  m{b)  <  fz{b)  <  M{b). 

Lemma  13  Under  condition  (Dl), 

(i)  there  exists  a  representation  Cjffj(i)  o/  tc/iere  the  r.v.’s  ^j, 

j  >  1,  are  independent,  identically  distributed 

(ii)  the  functions  gj{t),  j  >  1,  are  in  C\I); 

(Hi)  sup^tuh)axl\Dkiruiti,t2)  -  0  as  n oo,  0  <  *,  /  <  3. 

Proof,  (t)  The  Gaussian  process  Un{t)  —  Z)”-i  t  6  /,  is  centered  and  has 

covariance  function 


’■c/„(^i.f2)  =  Y,9j{h)9jit2),  h,t2  €  1. 
j=t 

For  each  tel,  the  sequence  {Unit} :  n  >  1}  is  Cauchy  in  L'^{Q,,A,P),  since,  for  aU 
p,  q  such  that  1  <  p  <  g,  we  have 

e(  E  «*(<))’  =  E  9j(<)  -  0, 

v=p+l  '  j=p+l 

asl<p<g-+oo.  Therefore,  this  sequence  converges  in  quadratic  mean  to  some 
Gaussian  r.v.  Ut  with  mean  zero  and  variance 

=  f^gjii)  =  ruit,t). 

j=l  jtzl 

Indeed,  - E^f^Ui{t)\  <  E^I\Ut-Un{t)f  ^  0  as  n  — >  cx).  Similarly,  for  any 

d  >  1  and  any  (ti, . . . ,  <<f)  €  P^,  the  sequence  {Unih), Un{td)  :  n  >  1}  converges 
in  quadratic  mean  to  a  Gaussian  R'^-valued  r.v.  Ut^,...,ti  ^^^h  mean  (0,. .  .,0)  and 
covariance  matrix  [ruiti,tk)]i<k,i<d- 
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The  family  of  the  distributions  of  the  r.v.'s  Ut^ . d  >  1, 

satisfies  the  symmetry  and  consistency  conditions  (Adler,  1981,  pp.  13-14).  By 
Kolmogorov’s  Theorem,  there  exists  a  Gaussian  process  U{t)  such  that  U{t)  -  Ut 
a.s.  for  all  t  €  I.  This  process  can  be  written  as  This  means  that  for 

each  t  ^  I,  U{t)  is  the  limit  of  quadratic  mean  (the  same  argument 

applies  to  (J7(ti), . . .  ,U{td.)))-  Moreover,  by  Cauchy-Schwarz  Inequality, 

\EiUih)U{t2))-E{Ur,{h)Un{i2))\ 

<  \E{U{h)U{t2))  -  E{Unitl)U(h))\  +  \EiUr,ih)U{h))  -  E(b\{h)Unit2))\ 

<  E^/^U\t2)E^/\Uih)  -  +  E^I'^Ul{h)E^I\U{t2)  -  Unih)?  -  0 

as  n  — »  00,  since  U[t)  is  the  limit  in  quadratic  mean  of  Un{i),  EU‘^{t2)  <  cc 
and  EUrXh)  -*  <  oo  as  n  ^  oo.  This  implies  that  E{U(ti)U(t2))  = 

lim„_oo  =  ru{ti,t2)- 

(ti)  From  (t),  we  deduce  that  gj{t)  =  E{U{t)^j)  for  all  j  >  1  and  all  t  6  /.  Besides, 
by  Theorem  2.2.2  in  Adler  (1981),  the  process  U{i),  i  £  I,is  4-times  differentiable 
in  quadratic  mean,  under  condition  (Dl).  Therefore, 

.  c-m)  i,) 

where  U{t)  denotes  the  quadratic  mean  derivative  of  U{t).  By  Cauchy-Schwarz 
Inequality  and  definition  of  we  have  gj(t)  =  E{U{t)^j).  Similarly,  we  can 
show  that  for  all  0  <  /:  <  4,  gj(t)  has  a  fe-th  order  derivative 

=  (24) 


where  denotes  the  fc-th  quadratic  mean  derivative  of  U{t).  Finally,  gj^\t) 

is  continuous  for  all  0  <  A  <  3. 

(iii)  For  fc  =  0,  since  gjit)  is  continuous  for  all  j  >  1  and  the  sequence  {Ej=i  9j{t)  ■ 
n>  1}  is  increasing,  it  follows  from  Dini’s  Theorem  that 


sup  3j(t)  =  snp 

<6/  j=n-H 


J=1 


0  (n  00 ). 


Since  (E|=n-(-i  <  'E,j=n+i9jih)'L%.n+i  9j{h)  for  all  g  >  n,  letting 

q  -*  oo,  we  obtain 


sup 


£  9jih)9j{t2} 


sup 


Mh,t2)  -^9j{h)9j{t2) 


0(n  -+  oo). 


For  Dkiru{ti,t2),  0  <  A:,  /  <  3,  by  Parseval’s  Inequality  applied  to  (24), 

OO 

J29f^\t)<EU^>‘^^t)<oo, 

j=i 
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for  all  f  €  /  and  all  0  <  fc,  /  <  4.  Therefore,  given  0  <  k,l  <  d,  we  can  construct 
a  Gaussian  process  ^  G  /,  as  in  (t).  Moreover,  since  the  functions 

J  >  1)  0  <  A:  <  3,  are  continuous  by  (ii),  by  Dini’s  Theorem, 

OO 

®  ^  °°)’  (25) 

for  0  <  A:  <  3,  and  the  covariance  function  of  this  process  is  the  uniform  limit  of 
as  n  — »■  00  for  €  /  x  /.  Finally,  it  results  from  (25)  that 

for  each  given  t2  €  /,  the  series  9jiti)9jih)  of  functions  of  ti  can  be  derived 
term-by-term  3-times.  Hence, 

OO 

Dkoru{h,t2)  -  Yl9f\h)9j{i2)-,  0<k<  'i.  (26) 

j=i 

Keeping  ii  constant  in  (26),  it  results  from  (25)  that  the  series  9j^\h)9j{t2) 
of  functions  of  t2  can  be  derived  term-by-term  3-times.  Hence, 

OO 

Dkiruih,t2)  =  0<  k,  I  <Z. 

Moreover,  by  Cauchy-Schwaxz  Inequality  and  (25)  the  series  !C>=i  fl'j*^(^x)5j^^(^2) 
converges  uniformly.  <> 

Remark.  Under  (Dl),  the  existence  of  functions  gj{t)  such  that  r[/(ti,t2)  = 
12%i  9.i(h)9j(t2}  results  from  Mercer’s  Theorem  (see,  e.g.  Adler,  1981,  Theorem 

3.3.1),  where  gj{t)  —  Xj  4>3{i),  j  >  1,  where  the  Aj’s  are  the  eigenvalues  and  the 
<^j(t)’s  are  the  eigenfunctions  of  the  integral  operator  with  kernel  ru(ti,t2). 

Let  Pn:  denote  the  orthogonal  projection  of  onto  IR”  x  {0, . . .},  defined 

by  P„x  =  (xi,...,in,0,...)  for  x  € 


Lemma  14  (i)  Under  condition  (Dl),  [|P„g(t)(|  =  1  -f  f„(t)  with  -*  0  as 

n  -+  OO  uniformly  for  t  £  I,  for  0  <  k  <3; 

(ii)  Under  conditions  (Dl)  —  (D2),  llPng'(<)ll  l|g^(0il  “uniformly  for  t  £  I  as 
n  OO  and  there  exists  Nx  such  that  inftg/  ||P„g'(t)|l  >  0  for  all  n  >  N\. 

Proof,  (j)  It  suffices  to  prove  a  similar  result  for  ||Png(0iP  =  J2%i9jW  =  1  ~ 
E^n+i  9]{t)-  For  each  0  <  A:  <  3,  the  A:-th  derivative  of  E^n+i  ffj  (0  is  a  linear 

combination  of  functions  of  the  form  9f^\^)9f^\^)  0  <  /i,l2  <  Ar.  By 

Lemma  13,  such  functions  converge  to  zero  uniformly  for  t  £  I. 

(ii)  The  first  statement  is  a  strciightforward  consequence  of  Lemma  13.  Since 
lig'(*)!l  is  a-  continuous  and  positive  function  of  t  £  I  and  I  is  compact,  6  = 
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inftg/ |ig'(f)||  >  0.  We  can  choose  Ni  such  that  sup^^/l  i|g'(OII  -  liPng'COII  I  <  ^/2 
for  all  n  >  Ni.  Then  infjg/  ||Png^(<)ll  >  ^/2  >  0  for  aJI  n  >  iVj.  0 

It  follows  from  Lemma  14  that 

P  s(t) 

gn{t)  =  t  6  I. 

is  well-defined  for  ail  n  >  N\  and  -*  uniformly  for  0  <  A:  < 

3.  The  functions  gn(0?  ^  parameterize  a  curve  7„  on  the  unit  sphere  of 

ffi”  X  {0, ...}  C  Moreover,  there  exists  N2  such  that  infig/ ||g^(t)||  >  0  for  all 
n  >  N2-  The  corresponding  unit  speed  parameterization  of  7„  is  defined  by  f„(s„)  = 
{/»i.i(«n),---,/n,n(5n),0,...)  =  g„(A“^(s„)),  where  s„  =  A„(t)  =  Jo  l|gn(*')ll<^^'  e 

Jn  =  [0,i„].  With  this  notation  Zn  =  supfg/A'„(t)  =  sup^^gj^  y„(s„)  where 
Xjt)  =  Tv^(t)Ur,(t),  UJt)  =  Ej=l  Ki(Sn)  =  KlSn)  " 

Ej=l  (jfnJ{Sn)  =  C^nCAnH-Sn))  and  Ty„(Sn)  =  Tx(XZ^{Sn)). 

Lemma  15  Under  conditions  (Dl)  —  (D2), 

(t)  An^^(t)  — >  A(*l(t)  uniformly  for  t  E  I,  for  0  <  fc  <  3; 

(it)  (A“^)(*')(s)  — >  (A~^)(*)(s)  uniformly  for  s  €  Jm  for  0  <  k  <3; 

(Hi)  fn‘\s)  -*  uniformly  for  s  £  Jn,  for  0  <  k  <3; 

(iu)  “*■  uniformly  for  s  €  Jn,  for  0  <  I  <3; 

(®)  c<,,„(s)  -»■  Cg(s)  uniformly  for  s  €  Jn- 

Proof,  (t)  By  Lemma  14,  g(i(t)  -*  g'(t)  uniformly.  Therefore,  A(,(t)  =  |(g'(0ll 
-^'(0  =  IlsXOil  uniformly  and  An(t)  =  fg  \\g\t')\\dt'  A(t)  =  ||g'(t'||dt'  uni¬ 

formly.  As  a  consequence,  =  An(T)  L  =  A(r),  Moreover,  the  fc-th  deriva¬ 
tive  of  A(,(t)  (respectively  A'(t))  is  the  product  of  a  polynomial  in  (gn‘^(t),gn*^(t)) 
(respectively  0  ^  hJ2  <  k,  and  ||g(,(t)||“^^^~^^  (respectively 

l|g'(t)||“^^^“^^).  Therefore,  the  result  follows  from  Lemma  14. 

(ti)  For  A:  =  0,  this  is  a  direct  consequence  of  (i).  For  A:  >  1,  the  result  follows 
also  from  (i)  since  the  A;-th  derivative  of  A“^(s)  (resp.  A"^(s))  is  an  expression 
involving  the  successive  derivatives  of  A„(t)  (resp.  A(t))  at  A~^(s)  (resp.  A~^(s)). 

(Hi)  This  is  a  direct  consequence  of  Lemma  13  and  (it). 

(iu)  This  follows  directly  from  (i),  (ii)  and  the  uniform  continuity  of  Tx(t)  and 
its  derivatives  over  the  compact  I. 

(u)  This  follows  from  (iii).  <) 

If  (Dl)  is  replaced  by  (D5)  in  Lemmas  14  and  15,  all  the  results  hold  for  0  < 
A  <  13  and  0  <  /  <  4. 
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Lemma  16  Under  conditions  (Dl)  —  (D4), 

(i)  There  exists  N3  such  that  for  all  n  >  N3,  6  =  ^ 

(ii)  there  exists  N4  such  that  for  all  n  >  N4,  the  curve  7„  parameterized  by  in(s)  ~ 
'^n('S)fn('S)  +  H,.(®)^n('S)>  ^  ^  hos  no  sclf-intersecUon. 

Proof,  (i)  The  proof  is  similar  to  the  proof  of  Lemma  14  (n). 

(ii)  We  first  show  that  there  exists  N4  and  €  >  0  such  that  fnC-s')  7^  fri(’S) 
for  all  $,s'  G  Jn,  [•s'  -  5|  <  e  and  n  >  N4.  By  a  quadratic  Taylor  expansion  of 
f„(s')  around  s,  f„(s  +  h)  =  +  hf^(s)  +  Knis,h),  where  ||R„(s, h)|| _<  Ch'^ 

and  C  =  sup„>jvj  sup^^j  ||f"(s)||  <  00  in  view  of  Lemmas  13-15.  Since  £^(s)  = 
+  ri;„(sj)f^(s)  +  T^Js)(Us)  +  as)),  lltoll  >  ry„(s)  +  r^Js)  >  ^  >  0  for 
all  s  €  J'n  whenever  n  >  N3.  Hence,  forn  >  N3,  ||f„(s4-h)-fn(5)ll  >  l^ll^-C'|h||  >  0 
for  aU  s  €  Jn  and  h  0,  jfti  <  t  =  6i{2C)  such  that  s  +  h  G  Jn- 

Let  Kt  =  {(^5 -s')  G  J  X  J  :  Is'  -  sj  >  e}.  Since  Kt  is  compact  and  the  function 
{s,s')  —*  ||f(s)  —  f(s')||  is  continuous  and  positive  on  K^,  64  =  inf(3^5')gj<-^  ||f(^)  " 
f(s')||  >  0.  By  Lemma  15  (iii)-(it)),  inf(,y)g/{-J|ii(s)  -  f„(s')j|  >  ^i/2  >  0  for 
n  >  N4.  The  conclusion  follows  for  n  >  N4  =  max{N4,N4).  0 

Lemma  17  Under  conditions  (Dl)  -  (D4), 

(i)  there  exists  Ns  such  that  for  all  n  >  Ns,  ^n,sis')  >  0  for  all  s'  ^  s; 

(ii)  there  exists  Ny  such  that  for  all  n  >  Nr  and  a.e.  s  G  Jn>  ^n(«).  fn(s0> 

fj^(s')  and  f^(s')  are  linearly  independent  for  all  s'  ^  s; 

(Hi)  there  exists  Ns  such  that  for  all  n  >  Ns  and  a.e.  s  G  Jn,  fnC'S),  f^(s),  fn(s'), 
^(s')  and  fni-s")  are  linearly  independent  for  all  s'  ^  s,  s"  ^  s,  s'  s". 

Proof,  (i)  By  a  Taylor  expansion  of  order  3,  it  follows  from  Lemma  15  (iii)-(iv) 
that  /3„,s(s  +  h)  =  (tv„(s)  +  T!fJ^s))h?/2  +  Rn(s,h),  where  Rn{s,h)  <  Cjlip  for 
some  constant  C  >  0.  Moreover,  by  Lemma  16  (i),  we  have  for  n  >  N'^  that 
PnAs  +  h)  >  h'^\6l2-C\h\\  >  0  for  all  s  G  Jn  and  h  7^  0  such  that  lft|  <  e  =  ^/(4C). 

Since  Ki  =  {(«,«')  G  J  x7  :  |s'-sj  >  c}  is  compact,  the  function  {s,s')  — »  ^s{s') 
is  continuous  and  positive  on  and  0n,s{s')  converges  uniformly  to  j8s(s')  for 
{s,s')  G  Jn  Jn  in  view  of  Lemma  15  (iii)-(it7),  it  follows  as  in  the  proof  of 
Lemma  16  (ii)  that  inf(,y)€Ke  i^n,»('S0  >  0  for  all  n  >  N'^.  The  conclusion  then 
follows  for  n  >  Ns  =  mBx(N^,  N^'). 

(ii)-(iii)  The  proof  of  (ii)-(iii)  is  similar  to  the  proof  of  Lemma  16  (ii).  First, 
using  a  Taylor  expansion  of  sufficient  order  and  Lemma  15,  we  show  that  the  Gram 
determinants  of  both  systems  are  positive  for  adl  s',  s'  ^  s  (resp.  (s',  s"),  s'  s, 
s"  ^  s  and  s'  7^  s")  sufficiently  close  whenever  n  is  large  enough.  Second,  we  take 
advantage  of  the  continuity  and  positivity  of  the  Gram  determinants  of  both  systems 
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of  sequences  for  |s'  -  s|  >  c  >  0  (resp.  |s'-  s|  >  e,  |s"  —  5|  >  f  and  |s'  -  5"|  >  e)  and 
of  the  uniform  convergence  of  the  corresponding  Gram  determinants  as  n  -t  oc.  0 

Lemma  18  Under  conditions  (Dl),  P{Z„  <  a)  -»  P{Z  <  a}  for  all  a. 

Proof.  First,  note  that 

IZ~Z„l<  sup  IX(i}  -  Xn(t)!  <  ^  .  sup  |r(t)  -  U„(i)l. 

tei  initg/rjv(Ej  «e; 

=  sup, ^!Var(U(t)  -  Un(t))  -  0  and  dl(ti,t2)  =  £((UiU)  -  f/„(t,))  -  (U(h)  - 
Un{l2))Y  ^  C\ti  —  t2p  for  n  large  enough  by  Lemma  13.  On  the  other  hand,  we 
have  the  following  inequality  due  to  Dmitrovskii  (Lifshits,  1986) 

P{sup \U{t)  -  i/n(<)|  >  ti}  =  2exp(-u'^/{2Kl))qniu), 
ta 

where  qn{u)  =  4.1  exp(2^^^6’P((/Cn/w)^'*^)(K:n/«)^^^(^(('«n/w)^^^))^^^  with  ^(Cn)  ~ 
€nlog^/^(Ci/€n)  when  €n  0  as  n  — >  oo  and  dn(ti,t2)  <  Clti  -  121-  Therefore, 
P{\Z  -  Z„|  >«}— »-0asn— >c»  and  converges  in  probability,  hence  weakly, 
to  Z.  Then  P{Zn  <  a}  — ♦  P{Z  <  a}  for  all  continuity  point  a.  As  P{Z  <  o}  is 
continuous  (Tsirel’son,  1975),  P{Zn  <  a}  — ►  P{Z  <  o}  for  all  o.  0 

By  Lebesgue’s  convergence  theorem  and  Lemmas  13-17,  it  follows  that  Mn{f>)  -* 
M(b)  for  all  b  and  m„(b)  m{b)  for  all  5  >  0.  Moreover,  M„(6)  <  M’'{b)  for  n 
large  enough  and  M{b)  <  M*{b)  where  M"{b)  =  Cq>{b0)  for  some  C  >  0  and 
^  >  0.  Hence,  the  sequence  {fz„}  is  weakly  compact  in  I^(R)  (Bourbaki,  1967,  pp. 
112-113),  which  means  that  there  exist  a  function  /  e  £^(I!t)  and  a  subsequence 
{/z„/}  such  that  fz^,{b)h{b)  db  -+  /(6)h(6)d6for  all  ft  €  X'”(R).  By  taking 

h{b)  =  l(_^,o](6),  it  results  that  <  fl}  — f{b)  db  for  all  a.  Therefore, 

by  Lemma  18,  Z  has  a  density  /z(6)  =  /(6).  Moreover,  the  same  result  holds  for 
any  accumulation  point  of  {/zn}-  Finally,  fz„{b)  -*  fz{b)  for  a.e.  b.  By  Tsirel’son 
(1975),  /z(5)  has  bounded  variation  on  every  compact  interval  of  R.  Hence,  fz{b) 
has  a  right  and  a  left  limit  at  each  point.  If  we  select  a  left  continuous  (say)  version 
of  it  follows  that  m(6)  <  fzib)  for  all  6  >  0  and  fz{b)  <  M{b)  for  all  6.  <> 
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A  Appendix 

Let  f(s)  =  (/i(5),...,/n(-5)),  s  &  J  -  [0,1-],  be  a  unit  speed  parameterization  of  a 
smooth  curve  7  embedded  in  the  unit  sphere  of  E”.  Therefore,  ||f(s)||  =  l|f'(>s)ll  =  U 
(f(s),f"(s))  =  -1,  ||f"(^)ll  >  0  for  all  s  €  J,  and  L  —  I7I  i:  'he  length  of  7. 
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At  each  point  M  =  M(s)  of  -r  we  can  define  the  unit  vector  tangent  to  -'i 
T  =  T{s)  =  f'(5)  and  the  principal  normal  vector  N  =  N(5)  =  f"(s)/|]r'(5)ji. 
Since  ||T||  =;  1,  N  =  T'/c  is  orthogonal  to  T  where  c  =  c(.s)  =  ||f"(s)||  defines  the 
curvature  of  7  at  the  point  M  =  M{s). 

For  each  s  £  J,  if  N(s)  ^  -f(s),  there  exists  only  one  unit  vector  K=K(s)  € 
Fect'^(f, T)  such  that  K  €  Vect(f,  N)  and  (N,K)  =  coso  >  0.  Moreover.  K(5) 
is  on  each  interval  on  which  (f,  N)  >  —1,  i.e.  c  ^  1.  If  c  =  1  we  can  define 
K(5)  =  K(s-).  Let  us  denote  by  Kj  =  Ki(s), . . K„_2  =  K„_2(s),  an  orthonor¬ 
mal  bcisis  of  V'ecf-‘-(f,T)  such  that  Ki  =  K  and  the  functions  K2,  • .  .,K„_2  are  ■ 
At  each  point  M  =  Af  (s)  of  7  the  moving  frame  (M,  f ,  T,  Ki, ....  Kn-2 )  is  orthonor¬ 
mal  and  the  matrix  of  (f', T',Kj, . .  .,K^_2)  with  respect  to  (f, T, Ki, - K„_2)  is 

antisymmetric: 

r  T  Ki  Ki  ...  k;_2 


where  D  is  a  (n  -  2)  x  (n  -  2)  antisymmetric  matrix,  Cg  =  Cg{$)  =  (Ki,T')  = 
c(Ki,N)  =  ccosa  >  0  defines  the  geodesic  curvature  of  7  at  the  point  M  =  M{s). 
By  definition  of  q,  we  have  (f,N)^  =  sin^a  =  1/c^  and  then  cj  -  -  1. 

Note  that  =  0  iff  cos  a  =  0,  i.e.  (f,N)  =  -1.  Finally,  it  can  be  shown  that 
K  =  (f  -f  n/Cg. 

If  g{t),  t  G  /  =  [0,r]  is  a  general  parameterization  of  7  with  ||g'(t)j|  >  0  for 
all  t  €  /,  the  corresponding  unit  speed  parameterization  of  7  is  given  by  f(.s)  = 
g(A“^(s)),  s  £  J,  where  s  =  A(t)  =  /J  ]|g'(u)|j  du,  t  £  I/i$  the  arc  length  of  7  from 
0  to  1.  We  have 
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